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Abstract 

Two identical ID autocatalytic systems with Gray-Scott kinetics — driven towards con- 
vectively unstable regimes and submitted to independent spatiotemporal Gaussian white 
noises — are coupled unidirectionally, but otherwise linearly. Numerical simulation then 
reveals that (even when perturbed by noise) the slave system replicates the convective 
patterns arising in the master one to a very high degree of precision, as indicated by 
several measures of synchronization. 

1 INTRODUCTION 

The beautiful talk in this Workshop by Jiirgcn Kurths [1] (as well as others dealing to some 
extent with the subject) exempts us to introduce the field of synchronization at large. Thus, 
hereafter we shall restrict our scope to the less explored subfield of the synchronization be- 
tween continuous systems [2-5] — concretely, to non-delayed synchronization between systems 
of stochastic partial differential equations. In particular, a topic that has been hardly addressed 
is the synchronization between noise-sustained structures (NSS) in systems undergoing a con- 
vective instability [6]. 

A convectively unstable regime is characterized by the fact that local perturbations to the 
steady state are advected more rapidly than their spreading rate [7] . When seen in a Lagrangian 
framework, the system is unstable; from an Eulcrian description, however, perturbations are 
"washed out by the flow". Macroscopic patterns named noise sustained structures (NSS) 
emerge in this regime if noise is present at all times. It is through dynamical amplification of 
random fluctuations that the system is driven out of its linearly unstable steady state towards 
the state sustaining NSS. Thus, if noise (or any external deterministic forcing) were not present, 
noncquilibrium structures could not arise. In fluid dynamics the NSS are a spatial macroscopic 
manifestation of amplified thermal fluctuations. 

NSS have been observed in fluid convection experiments (both in open flow configuration [8] 
and Taylor-Couette flows [9,10]), and their precursors have been also observed in nematic liquid 
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crystals [11]. They have also been numerically shown to exist in optical systems [6,12] (driven 
in this case by quantum noise) and recently, in a model autocatalytic chemical reaction — the 
Gray-Scott (GS) model — taking place in a differential-flow reactor [13]. 

The paper is organized as follows: In Sec. [2] a brief sketch is made of the GS model, and 
the rationale and features of the chosen master-slave coupling are pointed out. Section [3] 
introduces the details of the numerical integration scheme and discusses the features of the 
NSS arising in the uncoupled systems. Section [I] is devoted to a fairly thorough numerical 
characterization of the replication of NSS through complete synchronization. In particular, 
the behavior of the synchronization measures as functions of the parameters in the model is 
studied, and a numerical estimation is made of the robustness of the phenomenon. Finally, the 
main conclusions are summarized in Sec.[5l 

2 THE MODEL 

The GS model proposes three steps for the conversion of the precursor species P into the inert 
product C 
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The intermediate step has cubic autocatalytic kinetics. 

In the case we consider, the reaction takes place in a differential-flow reactor where A is 
immobilized, whereas B is free to diffuse and is also advected by the flow. Moreover, the 
reaction is maintained out of equilibrium by keeping the concentration of P constant (p = pa) 
and that of C zero (c = 0). Hence, the present model describes the irreversible decay of P 
towards a product C that is immediately removed from the reactor. 

After scaling concentrations by (A^/fci) 1 / 2 , time by k^ 1 and length by (-De/fe) 1 / 2 , the rate 
equations for system 1 read 

fi - a x b\ + y/ai~£i(r,t), 

V 2 6i-0V6i-6i+ai6?, (2.1) 

where /i stands for the scaled version of kapo and <f> for that of the fluid velocity v. The real 
Gaussian noise £i in the rate equation for a\ — with zero mean, variance o\ , and delta-correlated 
in space and time — accounts for fluctuations (either thermal or in p ). 

For n > 1, the uniform steady state (ai = = /i) becomes convectively unstable 

at some <fic{n), yielding to traveling waves with ±q c . Further details are found in [13] and 
references therein. 

Now we assume that system 1 drives another system (called thereafter system 2 and lying 
in a second differential-flow reactor) in a master-slave configuration. System 2 has the same 
values of n and <f) but its A-component is submitted to a spatiotemporal Gaussian white noise 
£2(1% t) with a possibly different noise variance a 2 : 

H - a 2 b\ + ■y/a 2 ^2(r,t), 

V 2 & 2 - <t> V6 2 - b 2 + a 2 bj + e (61 - b 2 ). (2.2) 
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e denotes the strength of the unidirectional linear coupling between both reactions. Besides 
being the simplest coupling that enables synchronization, it facilitates an approach to the 
stability analysis of the synchronization manifold [14]. 

3 NUMERICAL SIMULATION AND UNSYNCHRONIZED NSS 

We shall restrict hereafter to the ID case (the specificities found in higher spatial dimensions 



will be published elsewhere [14]). Equations (2.1l and (2.2) have been integrated using an 



Euler stochastic scheme in a grid of 16384 sample points with a grid space Ax = 0.1 and time 
step At = 0.0001. The parameters have been chosen as fi = 2.0, <fi = 9.5, o~\ — 10~ 7 . For each 
system, Dirichlct BC is assumed at the inlet of the reaction domain [a,i(0,t) = 6^(0, t) = 
/it, (i — 1, 2)] and Neumann BC at the outlet (x — L). The length L is chosen in such a way 
that spatiotemporal patterns develop well before they reach the outlet. 



For e = 0, Eqs. (2.1 ( — ( 2.2 1 describe two uncoupled reactions, identical with regard to the 
deterministic parameters but submitted to independent spatiotemporal noises which produce 
non-correlated NSS in both systems. These patterns have been characterized in Ref. [13]: they 
are dynamical structures that drift with the flow, disappearing on the right, whereas new wave 
excitations are continuously regenerated by dynamical amplification of noise. 

4 SYNCHRONIZATION OF NOISE-SUSTAINED STRUCTURES 

When e 0, some correlation is expected between the NSS in both systems. A handy measure 
of correlation for these snapshots are the deviation fields 

a(x, t) = ai(x,t) — a,2(x,t), 

f3(x,t) = b 1 (x,t)-b 2 (x,t). (4.1) 

Figures [TJi andfTj} are respectively snapshots of typical a\ and b\ profiles. The deviation fields 
are also depicted (in solid lines and in the same scales as a\, b\ respectively) in Figs.[l]i and[f|). 
The result is surprising, given that system 2 is also submitted to an independent spatiotemporal 
noise source: In the scales of Figs.[T] system 2 synchronizes completely to system 1. In other 
words, an effective replication of the NSS arising in the first reactor takes place at the second 
one, and time evolution — even under the influence of noise — does not spoil the high degree of 
synchronization. 

If we regard the system's evolution as a succession of snapshots like those of Figs, [I] natural 
quantifiers for this phenomenon (as functions of t) are the variances of a and (3: 



L 



a a (t) EE / [nV./l -{nrjtlr. 



(i 



ap{t) ee ^I^V(M)~</Wz, (4.2) 
with (if) ({) = j; J fix, t) dx (tp stands for a and j3 respectively), and the global synchronization 
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E ( t ) = \Tl [a 2 (x,t) + (3 2 (x,t)]dx. (4.3) 



Figure [2] is a plot of E, cr a and ap vs e for typical realizations (as stated before, the 
time evolution preserves the degree of synchrony). In the numerical simulation, (a) and 
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Fig. 1. A typical snapshot in the simulation of Eqs. (2.1 1 for e = 0.1. (a) ai vs x\ (b) bi vs x. 
Also shown — and in the same corresponding scales — are the deviation fields a and /3. The remaining 
parameters are n = 2.0, cj> = 9.5 and (T%,2 = 10~ r . 



remain below ~ 5 x 1CU 5 , i.e., E 2 ~ a 2 a + er| and basically E accumulates the information of 
both variances. On the other hand, correlations between the NSS remain during time evolution 
above .9999, indicating a very high degree of structure replication. In other words, the coupling 



in Eqs. (2.2) synchronizes the whole stochastic processes. 



A dependence of a a , ap (and thus of E) on the noise intensity 02 is to be expected. Figure[3] 
shows (here <j\ = 10 -7 ) that this is indeed the case, and maximum synchronization corresponds 
to (T2 = 0. In other words, only in system 1 does the noise play a constructive role (by pushing 
the system out of its unstable steady state). 

As usual, the synchronization between the stochastic fields a\ and ai (resp. h\ and 62) can 
also be viewed in the corresponding phase planes. As an illustration, Fig.[4]shows the dynamical 
correlation between b\ and 62 during the complete time history of a numerical simulation. 

One might wonder whether the proposed coupling is general enough, or whether the re- 
ported phenomenon is robust. To elucidate (at least partially) on these questions, we have per- 
formed numerical experiments where the coupling is switched on only after two well-developed 
and independent NSS are formed in each reactor. The results are shown in Fig. [5] The syn- 
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Fig. 2. Synchronization measures as functions of the coupling e: E (solid line), a a (dotted line) and 
073 (dashed line). The remaining parameters are /i = 2.0, (f> = 9.5 and 0-1,2 = 10~ 7 . 
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Fig. 3. E vs (12/01 for ai — 10 7 . The cr ajl a variances follow a similar behavior. The remaining 
parameters are n = 2.0, cj> = 9.5 and e = 0.1. 

chronization error decreases as soon as the coupling is switched on, a replication of system l's 
NSS takes place after a transient, regardless of the initial condition on reactor 2. In particular, 
there is no need to stabilize the second reactor prior to synchronizing it with the first one. 
This fact shows explicitly the robustness of the observed phenomenon, and suggests that the 
attraction basin of the synchronization manifold (02 = cti, 62 = b\) is large enough. 

In the same way, if the coupling is only switched on in a part of the reactors, numerical 
simulations (not shown) indicate that (even for 02 = 0) a replication of the "master structure" 
takes place in the coupled region of the slave system, opening the possibility of local replication 
(see below). 

We have also explored the synchronization of deterministic non-equilibrium structures (i.e, 
no noise in both systems). In particular, we have considered two cases of pattern generation: 



Prl-6 



s z - 




Fig. 4. (a): The (61,62) phase plane for the complete time history of a numerical realization of Eqs. 



( 2. 1 1 and (2.2 1 . A similar picture is obtained in the (01,02) plane. The parameters are = 2.0, 
6 = 9.5, 0-1.2 = 10~ 7 and e = 0.1. 




Fig. 5. E(t) when the coupling is switched on (to a value of e = 0.1) from t = 200 to t = 400 time 
units. The initial growth of E is associated with the pattern formation process in both reactors, which 
starts from the uniform solution perturbed by noise. Remaining parameters as in Fig. [2] 

(a) structures generated by the time evolution of an inhomogeneous initial condition (Fig. [HJ; 
and (b) structures generated by an oscillatory external forcing (Fig. [7]). In both cases, the 
numerical results (discused in the captions) clearly confirm the robustness of the proposed 
mechanism of synchronization. 



5 CONCLUSIONS 



By coupling unidirectionally (but otherwise linearly) corresponding points of two samples of 
the convectively unstable system under study, complete synchronization of macroscopic struc- 
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Fig. 6. Time evolution of a pulse-like initial condition for e = 0.1, without noise (<7i,2 = 0). The 
coupling is switched on between the vertical bold lines (1000 < x < 2000). Shown are the &i and /3 
fields for ti = 40, ti = 70, tz = 120, £4 = 160 and £5 = 250. Note that already for t±, the deterministic 
structure in system 1 is fully replicated in system 2, and they remain synchronized even without 
coupling. The other parameters are fj, = 2.0 and 4> — 9.5. 




Fig. 7. Synchronization of an externally forced deterministic structure. We show a snapshot of the 
61 and P fields. The left b.c. on system 1 [bi(x = 0,t) = y, + 0.5 sin(4.71t)] generates a traveling 
structure. The coupling (e = 0.2) is switched on at the right of the vertical bold line (300 < x) and 
the synchronization induced by coupling results evident. Remaining parameters as in Fig. [6] 

tures has been achieved, both for deterministic and stochastic dynamics. Figure [5] suggests 
that the synchronization attractor is very extended and a full replication of the structures is 
to be expected under very general conditions. This is a strong indication that the synchro- 
nization manifold (02 = ai, 62 — b\) is at least linearly stable. A complete stability analysis 
of the synchronization manifold will be published elsewhere [14]. We remark that the coupling 
synchronizes completely both systems (after a transient) regardless of the initial condition in 
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the "slave" system. Even more, the coupling may be defined in part of the system's extension. 
Thus, for (72 = 0, a (synchronized) structure arises in the slave system because of the coupling. 

We expect the phenomenon to have technological applications in the control of differential- 
flow chemical reactors, and eventually in the case where the convective structures in the master 
system are not noise-sustained ones, but carry useful information. 
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